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Abstract. We study the Einstein-Maxwell system of equations in spherically 
symmetric gravitational fields for static interior spacetimes. The condition for pressure 
isotropy is reduced to a recurrence equation with variable, rational coefficients. We 
demonstrate that this difference equation can be solved in general using mathematical 
induction. Consequently we can find an explicit exact solution to the Einstein- 
Maxwell field equations. The metric functions, energy density, pressure and the electric 
field intensity can be found explicitly. Our result contains models found previously 
including the neutron star model of Durgapal and Bannerji. By placing restrictions 
on parameters arising in the general series we show that the series terminate and 
there exist two linearly independent solutions. Consequently it is possible to find 
exact solutions in terms of elementary functions, namely polynomials and algebraic 
functions. 
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1. Introduction 

The first spacetime studied, modelling the interior of the relativistic sphere with uniform 
energy density, is the Schwarzschild interior solution which was found about ninety years 
ago. Since then considerable time and energy has been spent in finding exact solutions to 
the Einstein field equations for the interior spacetime that matches to the Schwarzschild 
exterior. The principal reason for this activity is that these solutions to the field 
equations for spherically symmetric gravitational fields are necessary in the description of 
compact objects in relativistic astrophysics. The models generated are used to describe 
relativistic spheres with strong gravitational fields as is the case in neutron stars. The 
detailed treatments of Stephani et al and Delgaty and Lake |2j for static, spherically 
symmetric models provide a comprehensive collection of interior spacetimes, satisfying 
a variety of criteria for physical admissability, that match to the Schwarzschild exterior 
spacetime. It is important to note that only a few of these solutions correspond to 
nonsingular gravitational potentials with a physically acceptable energy momentum 
tensor and a barotropic equation of state. A sample of the exact solutions to the field 
equations, which satisfy all the physical requirements for a relativistic star, is contained 
in the models of the Durgapal and Bannerji jHj, Durgapal and Fuloria Finch and 
Skea 5J, Ivanov [Bj, Lake [7j, Maharaj and Leach [H] and Sharma and Mukherjee 
amongst others. 

In the past most of the solutions found have been obtained in an ad hoc manner. 
We expect that a more systematic and formal study, such as the treatment of Maharaj 
and Chaisi ^U] for anisotropic matter, should lead to new classes of solution. Clearly 
there is a need to systematically study the mathematical properties and features of the 
underlying nonlinear differential equations. The analysis of John [H 3 indicates that 
reducing the condition of pressure isotropy to a recurrence relation with real, rational 
coefficients leads to new mathematical and physical insights in the Einstein equations for 
neutral matter. We attempt to perform similar analysis here in the coupled Einstein- 
Maxwell equations for charged matter. In this more general case we find that the 
condition of pressure isotropy leads to a new recurrence relation which can be solved in 
general. The Einstein-Maxwell system is important in the description of a relativistic 
star in the presence of an electromagnetic field. It is interesting to observe that, in the 
presence of charge, the gravitational collapse of a spherically symmetric distribution 
of matter to a point singularity may be avoided. In this situation the gravitational 
attraction is counterbalanced by the repulsive Coulombian force in addition to the 
pressure gradient. Consequently the Einstein-Maxwell system, for a charged star, has 
attracted considerable attention in various physical investigations including Mukherjee 
[12] and Sharma et al [T3] . 

In this paper we seek new exact solutions to the Einstein field equations, using 
a systematic series analysis, which may be used to describe the interior gravitational 
profile of a relativistic sphere. The approach produces a number of difference equations 
which we demonstrate can be solved from first principles. In section 2 we first express 
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the Einstein equations for neutral matter and the Einstein-Maxwell system for charged 
matter as equivalent sets of differential equations utilising a transformation due to 
Durgapal and Bannerji [B]. We choose particular forms for one of the gravitational 
potentials and the the electric field intensity, which we believe has not been studied 
before. This enables us to simplify the condition of pressure isotropy in section 3 to 
a second order linear equation in the remaining gravitational potential. We assume 
a series form for this function which yields a difference equation which we manage to 
solve using mathematical induction. It is then possible to exhibit a new exact solution 
to the Einstein-Maxwell field equations which can be written explicitly as shown in 
section 4. We consider two particular cases in section 5 which contain exact solutions 
found previously. In section 6 we demonstrate that it is possible to find two linearly 
independent solutions to the condition of pressure isotropy in terms of elementary 
functions by placing restrictions on parameters, that appear in the general solution, 
thereby permitting the series to terminate. In section 7 we express the general solution 
of the Einstein Maxwell system in terms of polynomials and algebraic functions. We 
briefly discuss some of the physical properties of our solutions in section 8. 



2. Field equations 

On physical grounds the gravitational field should be static and spherically symmetric 
for describing the internal structure of a dense compact relativistic sphere. Consequently 
we can find coordinates (t, r, 6, <fi) such that the line element is of the form 

ds 2 = _ e Mr) df 2 + e 2A(r) dr 2 + ^2 + ^2 Q &(j) 2y ^ 

For neutral perfect fluids the Einstein field equations can be expressed as follows 

4[r(l-e- 2A )r =p (2a) 



r 



' '1 - e' 2X ) + — e~ 2A = p (2b) 



- 2X 1 v" + v' 2 + - - v'X - - ) =p (2c) 




for the spherically symmetric line element ((TJ. The energy density p and the pressure 
p are measured relative to the comoving fluid 4-velocity u a = e~ u 5 a Q and primes denote 
differentiation with respect to the radial coordinate r. In the field equations (|2a )) -(|2~c f) 
we are using units where the coupling constant = 1 and the speed of light c = 1. 
The system of equations (|2aj) - (j2*c|) governs the behaviour of the gravitational field for a 
neutral perfect fluid. A different but equivalent form of the field equations is obtained 
if we define a new independent variable x, and new functions y and Z, as follows 

A 2 y 2 (x)=e 2v( - r \ Z(x) = e- 2X(r \ x = Cr 2 (3) 

so that the line element (P) becomes 

ds 2 = -A 2 y 2 dt 2 + -^-dx 2 + ^(d6 2 + sin 2 6d<p 2 ). 
4CxZ C 
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In the equations J3J) the quantities A and C are arbitrary constants. Under the 
transformation © the system ([2a j) -(|2"c j) has the form 

4*-* -§ <*•> 

4 Z i + ^± =L (44) 

y x C 

AZx 2 y + 2Zx 2 y + {Zx - Z + l)y = (4c) 

where the dots denotes differentiation with respect to the variable x. The set (J4~a|) -(|¥c j) 
is a system of three equations in the four unknowns p,p,y and Z. 

A generalisation of the system (|2~o j) - (j2"cj) is the Einstein-Maxwell field equations 
given by 

r(l-e- 2A )]' =P+ \ E2 (5a) 



-l(l-e-) + ^ e - = p-\ E 2 (56) 



2u' 1 

, . — & = P ~ -- 

- 2A ( v" + 1 y' 2 + -- v'X --) = p+-E 2 (5c) 
\ r r / 2 



1 



^e-\r 2 E)'. (bd) 



where E is the electric field intensity and a is the charge density. When the electric field 
E — then the Einstein-Maxwell equations (|5a j) -(|5d j ) reduce to the Einstein equations 
()2a|) - (f2~cj) for neutral matter. The system of equations ([5aj ) -([5*^ j) governs the behaviour 
of the gravitational field for a charged perfect fluid. If we use the transformation (J3j) 
then the Einstein-Maxwell system (I5ffl|) - fj5dj) becomes 

4*-* -h+rc 

( ■ E 2 x\ 
4Zx 2 y + 2Zx 2 y + ( Zx - Z + 1 - — J y = (6c) 

which may be easier to integrate in certain situations. 



3. Specifying Z and E 

We examine a particular form of the Einstein-Maxwell field equations ()6flj) -fl6 ^) by 
making explicit choices for the gravitational potential Z and the electric field intensity 
E. The system (I6"a j) - (j6dj) comprises four equations in six unknowns Z,y,p,p,E and 
a. By specifying the gravitational potential Z and electric field intensity E we are 
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in a position to integrate the condition of pressure isotropy (|6c|). The solution of the 
Einstein-Maxwell system then follow. We make the choice 

Z= 1 -±^ (7) 
1 + x 

where k is a real constant. In (JZj) we take k ^ 1. If = 1 then the metric function 
e 2A =l and the energy density is p = — 4p. To avoid negative energy densities, which 
are not physical for barotropic stars, we consequently take k 1. The choice 
was also made by Maharaj and Mkhwanazi Jl] and in their analysis of uncharged 
stars. Our objective is to confirm that this type of potential is also consistent with 
nonvanishing electromagnetic fields. Note that our choice contains, as a special case, 
the Durgapal and Bannerji |3| solution, which is widely applied as a relativistic model 
for neutral stars. Only the solutions for the cases k = | and k = — \ were documented 
previously for the uncharged case when E = 0. Other physically reasonable choices of 
the gravitational potential Z are possible; we have chosen the form (J2J) as it produces 
charged and uncharged solutions which are necessary for a realistic model. 
Upon substituting (JJJ) in equation (Jficj) we obtain 

/ E 2 (l + x) 2 \ 
4(1 + kx)(l + x)y + 2(k - l)y + I 1 - k K ! J y = 0. (8) 

As the differential equation (jHJ) is difficult to solve we first introduce the transformation 

1 + x = KX (9a) 

K = k -^ (96) 

Y(X) = y{x) (9c) 

so as to obtain a more convenient form. Substituting (|9"q |) -(f9"c j) in the differential 
equation (jHJ) we obtain 

x d 2 Y dY ( K 2 (l- K)E 2 X 2 \ , . 

AX(l-X)—-2- + (K + L—L. jy = (10) 

in terms of the new dependent and independent variables Y and X respectively. The 
differential equation (fTUj) may be integrated once the electric field E is specified. A 
variety of choices for E is possible but only a few are physically reasonable and generate 
solutions in closed form. We observe that the particular choice 



7712 _ — (~\~\) 

K 2 {l-K) X 2 1 J 

where a is a constant, simplifies (jl()[). The electric field defined in (fTT]) vanishes at the 
centre of the star, and remains continuous and bounded in the interior of the star for a 
wide range of values of the parameter k. Thus the choice for E is physically reasonable 
and it is a useful form to study the gravitational behaviour of charged stars. With the 
help of |TT| we find that |T0|) takes the simpler form 

«(i-x)g:-2§ + (* + Q )y = o. (12) 



aC KX -1 
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This is a special case of the hypergeometric equation. When a = the differential 
equation (fT^j) becomes 

«(l-X)g:-2§ + ^ = (13) 
and there is no charge. 

4. General series solution 

Since (J 12)) is the hypergeometric equation it is not possible to express the general solution 
in terms of elementary functions for all K + a. In general the solution will be given 
in terms of special functions. The representation of the solution in a simple form is 
necessary for a detailed physical analysis. Hence we attempt to obtain a general solution 
to the differential equation ()12p in series form. Later we show that it is possible to extract 
solutions in terms of algebraic functions and polynomials. 

Since X = is a regular singular point of the differential equation ()12|) we can apply 
the method of Frobenius about X = 0. We assume that the solution of the differential 
equation (|T2*|) is of the form 

oo 

Y=J2c n X n+r , c ^0 (14) 

71=0 

where c n are the coefficients of the series and r is a constant. For a legitimate solution 
we need to determine the coefficients c n as well as the parameter r. Substituting (fTlj) 
in the differential equation (|T2|) we obtain 

2c r[2(r-l)-l]X r - 1 + 

oo 

]T (2c n+1 (n + 1 + r)[2(n + r) - 1] - c n [A(n + r)(n + r - 1) - {K + a)]) X n+r = 0. (15) 

n=0 

The coefficients of the various powers of X must vanish. Equating the coefficient of 
X r ~ l in ()15j) to zero we obtain the indicial equation 

2c r[2(r - 1) - 1] = 0. 

Since cq ^ we must have r = or r = §. Equating the coefficient of X n+r in (fT5]l to 
zero we obtain 

A(n + r)(n + r-l)-(K + a) 

°n+i = -^r, — — r- Uo / — : — ; — 7i — c n, n > (16) 
2{n + 1 + r) [2(n + r) — 1\ 

which is the fundamental difference equation governing the structure of the solution. 

We can establish a general structure for the coefficients by considering the leading 

terms. The coefficients c±, C2, C3, . . .can all be written in terms of the leading coefficient 

c and we generate the expression 

_ ^ ^P + r)(p + r-l)-(K + a) 

n+1 M 2(p+l + r)[2(p + r)-l] °° [U) 

where the symbol fl denotes multiplication. It is possible to establish that the result 
()17|) holds for all nonnegative integers using the principle of mathematical induction. 
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We can now generate two linearly independent solutions, y\ and y%, from (|14j) and 
()17|) . For the parameter value r = we obtain the first solution 



Y 1 = c 



4p(p 



(K + a) 



55 2(p+l)(2p-l 



n+1 



4p(p- 1) - (K + a) + x 



1+ £S> 2(p+l)(2p-l) 



K 



n+1 



For the parameter value r = | we obtain the second solution 



Y 2 = c X 



1 + g TT (2p + 3)(2p+l)-(K + a) xn+1 



n=0p=0 



(2p+5)(2p + 2) 



Vi = c 



1 + X\ 2 



~ » (2p + 3)(2p + 1) - (X + a) f 1 + ^ " +J 



(2p + 5)(2p + 2) 



(19) 



for the choice pijl . is given 

(20) 



n=0p=0 

Thus the general solution to the differential equation 
by 

y = ayx(x) + by 2 (x) 

where a and b are arbitrary constants, K = and y\ and y 2 are given by (|TH|) and 
(|19|) respectively. By inspection it is clear that y\ and y% are linearly independent 
functions. From (p3aj) - (l6dj) . (fTSj) and (fT9*j) the general solution to the Einstein-Maxwell 
system becomes 

1 + x 



e 2A 

P_ 

C 

P_ 

c 

El 
c 



1 + kx 
A 2 y 2 

(l-fc)(3 + x) 

(i + xY 



akx 



2 1 



a foe 



(21a) 
(216) 
(21c) 

{21 d) 

(21e) 

where y = ayi(x) + 6y 2 (a;). We believe that (|21aj) - (l21ej) is a new solution to the Einstein- 
Maxwell field equations. 



^l + fcap j/ | (fc-1) | 



(1 + x) y 
afor 



1 + 2(1 + x) 2 



5. Particular cases 



From the Einstein-Maxwell solution (j21aj) - ()21dj) we can generate a number of physically 
reasonable charged and uncharged models for particular choices of k and a. If we set 
a = then 

1 + x 



,2A 



„2i/ 



1 + foe 

AY 



(22 a) 
(226) 
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p_ (1 -k)(3 + x) 

C ~ (1+x) 2 

V_ = 1 (1 + kx) y (k - 1) 

C (1 + x) y (l + x) 



(22 c) 
(22 d) 



which corresponds to a neutral relativistic star. We believe that the uncharged solution 
()22aj) - fl22flf is also a new solution to the Einstein field equations (|4*fflj) - fj4*cj) . It does not 
appear in the comprehensive list of solutions presented by Delgaty and Lake 2j. In 
the solutions ()21a}) - (l21e|) and ()22aj) - fl22dj) the gravitational potentials A and v are well 
behaved. Clearly the energy density p is positive at the origin if we choose k < 1. The 
pressure is finite at the origin. These are desirable features in a stellar model. 

When K + a = 3 the series in (fTH|) terminates. It is then possible to write the exact 
solution to the Einstein-Maxwell system in terms of elementary functions. The explicit 
form of the solution is given by 
(K — 1)(1 + x) 



-,2A 



C 

P_ 
C 



(K-l-x) 
A 2 [ci(l + x)3 +c 2 (K-l 
2K(3 + x) + (3 - K)x 
2(K-l)(l+x) 2 



x) 



■■\K + 2 + 2x) 



(23 a) 
(236) 
(23 c) 



K - 1 

Ci(l + x) 



[5K — 6 — (K + 6)x] + c 2 {K - 1 - x) 5 [4K - 12 - K 2 - (12 + 2K)x] 



ci(l + x)i + c 2 (l + x)(K - 1 - x)*(K + 2 + 2x) 



+ 



(3 - K)x 



El 
~C 



2(1 ~K)(1 
(3 - K)x 



x) 



(23d) 

[l-K)(l + x) 2 (23g) 
where K = (k — l)/k and x = Cr 2 . Clearly (j23aj) - ()23e|) is a special case of the 
general solution ()21aj) - fl21e|) . The solution (j23aj) - (l23ej) has the advantage of being given 
completely in terms of elementary functions which makes an analysis of the physical 
features of the model possible. When K = 3 (i.e. k = — |) and a = we obtain 

2(l+x) 



-,2A 



'2;/ 



(2-x) 



A 2 



ci(l + x)* + c 2 (2 - x)a(5 + 2x) 



P_ 

C 

P_ 

c 



3(3 + x) 
2(1 + x) 2 



ci(l + x)2(l -x)- c 2 (2 - x)*{l + 2x) 



(24a) 
(246) 
(24c) 

(24a0 



ci(l + as) 5 + c 2 (l + x)(2 - x)a(5 + 2x) 

for an uncharged relativistic stellar model. The special case (I24aj) - fl24fl| is the same as 
the result of Durgapal and Bannerji j2] and Maharaj and Mhkwanazi [Tl]. We point 
out that Maharaj and Mhkwanazi ^3] had a numerical mistake in their calculation of 
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the pressure p which has been corrected in our solution. We believe that the solution 
()23aj) - (j23dj) is important in the study of charged stars as it contains the Durgapal and 
Bannerji jHj model which has been shown to be consistent with a realistic dense star. 
Extensive studies of the Durgapal and Bannerji solution, as indicated in the compendium 
by Delgaty and Lake has proved that all the criteria for physical acceptability are 
satisfied. It is consequently used in many astrophysical studies that model neutron 
stars. 



6. Terminating series 

The general solution (|2T)|) can be expressed in terms of the special functions, namely 
hypergeometric functions. For particular values of K and a the series solution can be 
given in terms of elementary functions as demonstrated in section 5. This is possible in 
general because the series (|T%|) and (fTTIJ) terminate for restricted values of the parameters 
K and a. Using this feature we obtain two sets of general solutions in terms of 
elementary functions, by determining the specific restriction on K + a for a terminating 
series, as demonstrated in the following sections. 



6.1. The first solution 

On substituting r = in equation (|TKJ) we obtain 



Ci+l 



4i(i -l)-(K + a) 



Ci, i > 0. 



(25) 



(2i + 2)(2i-l) 

If we set K + a = An(n — 1), where n is a fixed integer, then c n+ i = 0. Clearly the 
subsequent coefficients c n+2 , c n+3 , c n+4 , . . . vanish and equation (j23|) has the solution 



q = 4n(n — 1) 



(-A) l - 1 (2i - l)(n + i - 2)! 



{2%)\{n-i)\ 

Then from the equations (fTfj) (when r = 0) and 



c , 1 < i < n. 



we obtain 



. . M,-,)X ( " ir ■ {2/ : [){ ": / 2)! .v 



i=l 



(2i)\(n 



iV. 



(26) 



(27) 



where K + a = An(n — 1). 



On substituting r = | in (fTB^) we obtain 



(2z + 3)(2z + 1) -(K + a) 



Ci, i > 0. 



(2z + 5)(2z + 2) ~" (28) 

If we set K + a = (2n + 3)(2n + 1), where n is a fixed integer then c n+ \ = 0. Also 
the subsequent coefficients c„+2, c n+ 3, c n+ 4, . . . vanish and equation (|2~£j) can be solved 
to yield 

3(-4) l (2i + 2)(n + i + 1)! 



(n + l)(n-i)!(2i + 3)! 



-c , 1 < i < n. 



(29) 
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Then from the equations (J 14)) (when r = ~ ) and (|29|) we obtain 



Yt = coXi 



3 " (-4)'(2» + 2)(n + » + !)! ■ 
{n + l)h (n-i)\(2i + 3)\ 



(30) 



where K + a = (2n + 3)(2n + 1). The polynomial and algebraic functions (|27j) and 
(|3U|) comprise the first solution of the differential equation ((TJ)) for appropriate values 
of K + a. 

6.2. The second solution 

We can use the form of the particular solution in section 5, expressed in terms of 
elementary functions, to simplify the representation of the second solution. The special 
solution in section 5 contains terms of the form (1 — X)^(l + 2X) which is product of 
(1 — X)z and a polynomial. This suggests that the second solution in general is of the 
form 

Y 2 = (l-X)*u(X) 

where u(X) is an arbitrary function. We now take Y 2 to be the generic second solution 
of (JT2J) and explicitly determine u(X). Equation (fT^j) gives 

4X(1 - X)il - 2(1 + 2X)u + (1 + K + a)u = (31) 

where dots denote differentiation with respect to X. 

Since X = is a regular singular point of the differential equation (}3~Tj) we can 
apply the method of Frobenius. We assume that the solution is of the form 

oo 

u=Y,CnX n+ \ C ^0. (32) 

n=0 

On substituting ()32|) in the differential equation (|31|) we obtain 
2c r[2(r - 1) - l]X r - 1 + 

oo 

]T (2c n+1 (n + 1 + r)[2(n + r) - 1] - c n [A(n + rf - (1 + K + a)]) X n+r = 0. (33) 

n=0 

The coefficients of the various powers of X have to vanish. Setting the coefficient of 
X r ~ 1 in ()33|) to zero we obtain the indicial equation 

2c r[2(r - 1) - 1] = 0. 

Since Cq ^ we must have r = or r = | as in section 4. Equating the coefficient of 
X n+r in to zero we obtain 

= 4(n + r)'-(l + g + a) 
+ 2(n + r + l)[2(n + r)-l] V ; 

which is the basic difference equation governing the structure of the solution. 

We establish a general structure for the coefficients by considering the leading terms. 
On substituting r = in equation (|3^jl we obtain 

4i 2 -(l + K + a) 

Q+1 " (2t + 2)(2t-l) Q ' { } 
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We assume that K + at = (2n + 3)(2n + l) where n is a fixed integer. Then c n+ 2 = from 
(f33|) . Consequently the remaining coefficients c n+3 , c n+4 , c n+5 , . . . vanish and equation 
(jHSJ) has the solution 



d = 4(n + r 



(-4) i - 1 (2i - l)(n + i)\ 



Co, 1 < i < n + 1. 



(2z)!(n-z + l)! 
Then from equation ()32j) fwhen r = 0) and (JBTIj) we obtain 

L , M , 1 ^ 1 (-4y- 1 (2»-l)(n + »)! y< 
u = Cq 1 + 4(n + 1) / , : — A 

i=l 

Hence we have the result 



(36) 



(2i)!(n 



1)! 



r 2 = Co (i-x)- 



"i + 4(, + i)E ( - 4fl(2 ^ 1)(r2 + 2)! 

i=l 

where if + a = (2n + 3)(2n + 1). 



On substituting r = | in equation 



(2i)!(n-i + l)! 
we obtain 



X 1 



(37) 



(2z + 3) 2 -(l + if + a) 



-a. 



(2* + 5)(2* + 2) ^ (38) 
We assume that K + a = 4n(n — 1) where n is a fixed integer. Then c n _i = from ([38)1 . 
Consequently the remaining coefficients c n , c n+ i, c n+2 , • • • vanish and (J3~%j) can be solved 
to yield 

3(-4) i (2i + 2)(n + i)! 



n(n-l)(2i + 3)!(n-i-2)! 



c , i < n 



(39) 



Then from the equations (|32|1 (when r = |) and ()39|1 we obtain 

-4)*(2i + 2)(n + i)! 



Hence we have the result 



1 + 



n-2 



T)5 



r 2 = c (i-x)2X 



;/(;/ - l) fr{ (2i + 3)!(n-i-2)! 
3 



1 + 



g (-4) J (2z + 2)(n + i)! 



(40) 



n(n-l) f^ x (2i + 3)!(n-i-2)! 
where if + a = 4n(n — 1). 

The solutions (|3*Tj) and (jjQj) generate the second solution of the differential equation 
(fT2Jl which are clearly independent from the solutions (}2T|) and (j3~Uj) . The quantities (}3*T|) 
and (|4())l are products of polynomials and algebraic functions. 



7. Elementary functions 

Thus we have generated general solutions to the differential equation (fT2*|) by restricting 
the values of K + a so that polynomials and product of polynomials with algebraic 
functions are possible as solutions. Collecting these results we have the first category of 
solutions 

- 1+4(n+1 )g<- 4 );::' 2 '- i ""i') ! x- 



Y = a(l-X)- 



i=i 



(2i)!(n-i + l)! 



1+ 



n + 



" (-4)<(2» + 2)(n + » + !)! , 
l)h (n-i)\(2i + 3)\ 



(41) 
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for K + a = (2n + 3)(2n + 1) where a and b are arbitrary constants. In terms of x the 
solution (j4*Tj) becomes 



y 



{K-l-x 
a I 

+ b 



K 

1 + x 
K 



1 + 



1 + 4(n 
3 



(_ 4 )- 1 (2z- l)(n + i)! fl + z 

i=i 



E 



(2i)!(n-i + l)! 
-4)*(2i + 2)(n + i + l)! 



+ (n-i)!(2i + 3)! 

The second category of solutions is given by 

3 



1 + x 
K 



(42) 



y = a(l-X)2j 



1 + 



n[n 



n-2 

E 



4) l (2i + 2)(n + z) 



+ b 



1 + 4n(n - 1)5Z 



-1) ^ (2z + 3)!(n-z-2 
-4) i ~ 1 (2i - l)(n + i- 2) 



(43) 



(2i)!(n-<)! 

for K + a = 4n(n — 1) where a and 6 are arbitrary constants. In terms of x the solution 
becomes 



K-l-x 



+ b 



K 

1 + 4n(n 



1 + x 



K 

n 

i)E 



i + 



n(n — 1) 
(_4) i - 1 (2i - l)(n 



E 



(-4)*(2i + 2)(n + i)! + s 



L (2i + 3)!(n-i-2)! 
z - 2)! /I + x s 



7\ 



0« 



(44) 



(2i)!(n 

It is remarkable that the solutions (|4*2*|) and (|44jl are expressed completely as 
combinations of polynomial and algebraic functions. It is rare to find general solutions, 
considering the nonlinearity of the gravitational interactions, to the field equations in 
terms of elementary functions. We have expressed our solutions in the simplest possible 
form. This has the advantage of simplifying the analysis of the physical properties 
of the dense star. Observe that our treatment has combined both the charged and 
uncharged cases for a relativistic star. If we substitute a = in and then we 
can obtain the solutions for the uncharged case directly. Consequently our approach 
has the unexpected but very desirable feature of producing an uncharged (possibly new) 
solution to equations (jjaj)-(jli| from the charged solutions when E = 0. We believe that 
the solutions obtained in this paper to the Einstein ()4a|) - fl4*cj) and Einstein Maxwell 
(jf)ojl - f)6djl field equations have not been found before. 

From our general class of solutions (J4*2*|) and (jUj) it is possible to generate particular 
solutions found previously. If we take K = 3 and a = (n = 0) then it is easy to verify 
that the equation (J4*2~j) becomes 



y = Cl (2 - 2)2 (5 + 2x) + c 2 (l + X 



where C\ = a/9 and C2 = b/3? are new arbitrary constants. Thus we have regained the 
Durgapal and Bannerji j3J neutron star model. Other explicit functional forms for y 
are obtainable which could be useful in applications for a dense star. As an example 
suppose that K + a = 8(n = 2) then from (jUj) we obtain 

y = Cl (K - 1 - x)3(l + a;)i +c 2 \K 2 + AK(1 + x) - 8(1 + x) 
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where C\ = a/K 2 and c 2 = b/K 2 are new arbitrary constants. It is now possible 

to generate an exact solution to the Einstein-Maxwell system (|6o jl -(|6d jl in terms of 
elementary functions when K + a = 8. This is given by 

e- = (gj^Kllg) (45a) 

e 2u = A 2 [c x {K - 1 - + x )l + c 2 (fsT 2 + 4K"(1 + x) - 8(1 + a;) 2 )] 2 (456) 
p 2K(3 + x) + (K-8)x 



C 2(K-l)(l+x) 2 
C 



(45 c) 



2{K -l-x)^ 


Ci(l + 2^ 


5(3K- 


4(1 + x)) + 8c 2 (K - 1 - x)%(K - 4(1 + x))" 


(K-l)(l + x) 




1 — x)a 


(1 + x)§ + c 2 (iT 2 + 4A"(1 + x)- 8(1 + x) 2 )" 



(K-l)(l + x) 2(1 -K)(l + x) 2 
E 2 (K - 8)s 



(45 e) 



C (l-K)(l+x) 2 

where K — (k — l)/k and x = Cr 2 . The solution (j45a|) - (l45ej) is given in a simple form 
which facilitates a physical analysis. 

8. Discussion 

We have found new solutions ()21aj) - (j21ej) to the Einstein-Maxwell system using a 
systematic series analysis that produces a number of difference equations which can 
be solved in general. A useful feature of the approach is that we regain the Durgapal 
and Bannerji neutron star model j2j as a special case which suggests that our class of 
solutions are physically reasonable. We briefly consider some physical features of the 
solutions of interest. 

Firstly, in the general solution Ij21qj) - (l21e|) . when studying models of charges 
spheres, we should consider only those values of k for which the energy density p, 
the pressure p and the electric field intensity E are positive. The choice of k must 
ensure that the gravitational potential e 2A remains positive; the remaining potential 
e 2v is necessarily positive. Clearly a wide range of charged spheres, with nonsingular 
potentials and matter variables, are possible for relevant choices of k. The interior 
metric ((TJ must match to the Reissner-Nordstrom exterior spacetime 

, , f 2M Q 2 \ 2 / 2M Q 2 Y\ o 9/in9 
ds 2 = - 1 + -V dt 2 + 1 + ^t) dr 2 + r 2 (d9 2 + sin 



across the boundary r = R. This yields the relationships 

~a yi (CR 2 ) + by 2 (CR 2 ) 



2M , Q 2 __ a2 




1 + CR 2 
1 + kCR 2 
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between the constants a, b, k, A and C. This shows that continuity of the metric 
coefficients across the boundary of the star is easily achieved. The matching conditions 
at the boundary may place restrictions on the function v and its first derivative (for 
uncharged matter) and the pressure may be nonzero (if there is a surface layer of charge); 
there are sufficient free parameters available to satisfy the necessary conditions that may 
arise from a particular physical model under consideration. 

Secondly, we observe that our solutions may be interpreted as models for anisotropic 
spheres (which may be charged or uncharged) where the parameter a plays the role of the 
anisotropy factor. The solutions found depend smoothly on the parameter a; isotropic 
and uncharged solutions can be regained for a = 0. For recent analyses of the physics of 
anisotropic matter see Chaisi and Maharaj ^2], Dev and Gleiser [TJi], and Maharaj 
and Chaisi [TO] . 
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